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1 Introduction 

The forthcoming luminosity upgrade of the Fermilab TEVATRON pp collider up to least 500 pb -1 
per year and the installation of the precise tracking systems in the existing detectors allowing the 
detection of the secondary vertices from b-quarks ( b-tagging) will open new opportunities for the 
Higgs boson search at the TEVATRON. 

In [1, 2] it was shown that at the TEVATRON energies (y/s = 2 or 4 TeV) the most promising 
will be Higgs boson search in the processes of Higgs production in association with the electroweak 
W and Z bosons. These processes are the direct analogue of the bremsstrahlung mechanism of the 
Higgs production in e + e~ collisions. Calculations of the Higgs signal reactions 



have been done for LHC and SSC energy range in a number of papers [3]. In the paper [1] the 
calculation of the reaction (1) as well as complete tree level calculation of the reaction pp — > Wbb-\- X 
have been done at TEVATRON energies. In the paper [2] the reactions (1) and (2) at TEVATRON 
energies have been studied including QCD radiative corrections. 

In this paper we present the more detailed study of the Higgs signal processes in the Higgs mass 
range Mjj < 140 GeV where the Higgs decay into a pair of bb quarks dominates, as well as the study 
of the basic background processes including an analysis of the 4-fermion final state distributions. 

The paper is organized as follows. In part 2 we will present the results for the total rate of 
the 4-fermion final state reactions. The rates, most promising signatures and various kinematical 
distributions of the Higgs production processes (1,2) are analyzed. In part 3 we present the results of 
the background simulation. The optimal kinematical cuts are proposed based on the comparison of 
the signal and background distributions. The effective mass distribution of the bb-pair after applying 
cuts and the extraction of the Higgs signal from the background are discussed. 

2 Total event rate and the Higgs boson contribution 

Due to the fact that the dominating decay mode of the Higgs in the mass interval Mjj < 140 GeV 
is H — > bb the Higgs production reactions (1), (2) lead to the following parton final states: 



pp -> W ± H + X, 




pp 



ZH + X 




pp 



pp 



pp 



bbfiv + X 
bbvv + X 
bbl+r + X 



(3) 
(4) 
(5) 
(6) 
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where / is one of the charged lepton e, //, r. 

In the paper we will consider in detail the reactions (3)- (5) which correspond to the leptonic 
decay modes of W- and Z-bosons. These reactions provide more clean event signatures for detection, 
namely 



pp -> 2b-jets + l+ + X (7) 

pp ^ 2b-jets+ Mfi±) + X ( 8 ) 
pp -> 2b- jets + /+/"+ X, (9) 

where j)±(/E±) denotes missing transverse momentum (energy). 
The reaction (6) leads to the signature 

pp -> 2b- jets + 2- jets + X (10) 

with the largest event rate. However the reaction (6) is more complicated phenomenologically be- 
cause the complete 2 — > 4 calculation has to be done in this case, and experimentally because the 
corresponding signature (10) is more difficult for separation of the Higgs signal contribution from the 
background one. A detailed study of the reaction (9) will be done in a separate paper. 

The complete set of tree level Standard Model diagrams contributing to the final states (3), (4), 
(5) includes the following numbers of 2 — > 4 diagrams: 

18 diagrams for the parton reaction qq — > bbl^v; 

33 diagrams for the reaction qq — > bbl + l~; 

16 diagrams for the reaction gg — > bbl + l~; 

15 diagrams for the reaction qq — > bbvv\ 

8 diagrams for the reaction gg — > bbvv , 
where / = e, //, r and q = u,d,s quarks. 

We will describe a calculation procedure in the next part but here we will formulate the results. 
The total cross section behavior is presented in Fig. 1 at yfs = 2 and 4 TeV versus the Higgs boson 
mass. For comparison the contribution of the Higgs signal diagrams ( processes (1), (2) ) is presented 
also in the same Fig. 1. One can see that, for instance, at y/s = 2 TeV the Higgs production cross 
section is about 0.3 - 0.03 pb depending on the Higgs mass and therefore one can expect 300 - 30 
signal events for integrated luminosity L = 1000 pb -1 . However the total reaction rates are about 
two orders of magnitude larger that the Higgs signal. That's why one has to make a detailed analysis 
of different distributions to find out whether it is possible to apply some set of cuts to suppress the 
huge background. 

Higgs production event rates with the final states (3)- (5) are shown in Table 1 for the ener- 
gies y/s = 2 TeV and y/s = 4 TeV. The number of events corresponds to the integrated lumi- 
nosity L = 1000 pb -1 . The results for the signal event rates have been obtained by means of 
PYTHIA 5.6/JETSET 7.3 Monte-Carlo generator [4] and include QCD radiative corrections and the 
effects of initial and final state radiation. 

Various H° signal event distributions for the both mechanisms (1), (2) are presented in Fig. 2 for 
the case M H = 80 GeV and y/7s = 2 TeV and in Fig. 3 for M H = 120 GeV and y/7s = 4 TeV. In 
particular we consider the b-quark and charged lepton transverse momentum distribution — p±b and 
P±h the total bb pair transverse momentum distribution — P]_bli an d missing transverse momentum 
distribution — fp^. The following detailed discussion in section 3 leads us to the conclusion that 
application of the appropriate cuts on these variables should suppress the overwhelming background 
to a level small enough. 
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It is important to point out, that in some cases it is not possible to distinguish signature (7) 
connected with WH process (1) and signature (8), (9) connected with ZH production mechanism 
(2). Indeed in the case if a lepton stays in the beam pipe or missed in detector, signature (7) looks 
like (8) one. If the only one lepton is detected, the signature (9) imitates signature (7). Due to this 
fact it is reasonable to consider the Higgs signal contribution together. Of course, in this case one 
needs to simulate the corresponding background as well. 



3 Background study. Extraction of the Higgs signal. 

From the complete set of diagrams mentioned above the dominant contributions to the reactions 
(3)-(5) come from several parton subprocesses. In Fig. 4. the diagrams for QCD 2 — > 3 type main 
background subprocesses 



qq — >■ Wbb 
qq — > Zbb 

99 ^ Zbb (11) 
are presented. 

A very important background is connected with 2 — > 2 electroweak subprocesses 



qq -»■ WZ 

qq - (Z/7*)(Z/ 7 *) (12) 
with subsequent leptonic decays of W- and Z-bosons or virtual 7*'s. 

An important contribution to the background comes from the parton reaction of single top quark 
production 

qq — > tb or tb (13) 

with subsequent decay of top quark t — > Wb. In our calculations we have put the top mass equal to 
170 GeV. The contribution of the background (13) has not been taken into account in the previous 
studies [1, 2]. Below we will see that this contribution becomes more important with the growth of 
the invariant 66-pair mass. 

Contributions of the other diagrams as well as different interference terms from the complete set 
of Feynman diagrams are small ( about 3-5% ) and one can neglect it in the analysis. The situation 
here is very similar with 4-fermion final states reactions in e + e~ collisions e + e~ — > bbfi + fi~ } bbe + e~ } 
bbvv [5], where also the main contribution comes either from resonant diagrams or from diagrams 
with soft photons. The reactions (11) were considered in the literature ( Wbb in [6], Zbb in [7]) as 
well as processes (12) which were calculated with the next to leading order radiative corrections ( 
ZZ in [8] and WZ in [9]) giving an effect of order 30%. 

We present the results of the calculations taking into account decays of W- and Z- bosons with 
four body final state kinematics. The following procedure has been used in the case of 2 — > 3 sub- 
processes (11). Feynman diagrams generation, analytical calculation of the squared matrix elements 
and generation of the corresponding FORTRAN codes have been done with the help of CompHEP 
software package [10]. Based on FORTRAN code, a Monte-Carlo generator for event flow at the 
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parton level was made. This generator has been written as a user process subroutine for PYTHIA 
package which develops the subsequent decay and fragmentation of partons. One of the processes of 
our generator gg — > Zbb was already in PYTHIA, and it has been used as an additional test. For 
the subprocesses of the 2 — > 2 type (12) the standard PYTHIA built-in processes have been used for 
event simulation. 

The HMRS set from PDFLIB [11] for the proton structure function have been chosen. However in 
order to study the level of uncertainties of calculations due to different parameterizations two other 
structure functions (M-T set Bl and GRV LO) have been used also. As a result we have found that 
the level of uncertainty is about 25 — 30%. 

In Fig. 5,6 the results of calculation and simulation of the same distributions as for the signal are 
presented for background processes: 

pp — > WZ — > lvhh\ 

pp — > Wbb — > Ivbb] 

pp — > ZZ — > bbll^p); 

pp — > Zbb — > bbll^p) 

at y/s = 2 TeV and y/s = 4 TeV respectively. The main differences in the distributions are connected 
to the fact that the events of QCD backgrounds (11) (Fig. 4) are concentrated most of all in the region 
of small p±. As a result of the detailed comparison of the signal and backgrounds distributions the 
appropriate set of the optimum cuts has been found: 

• p ±b > 20 GeV; p n > 15 GeV ; 

. p Lbb >20GeV ; (14) 
the following cuts have been used also for the event selection typical for TEVATRON detectors [2]: 

• jp L > 20GeV 

• 1 7/ 6 1 < 1.5; \r]i\ < 2 



. |Ai^| > 0.7; \AR bl \ > 0.7 (|Ai^| = y/A^ + 

• resolution for M bb : AM b - b /M b - b = 0.8/ \/m^ +0.3, 

which based on jet energy resolution 

^ 2AM bb + M bb > M bb > -2AM bb + M bb . (15) 

This set of cuts reduce the total background more than 40 times. We would like to stress the efficiency 
of p^i cut which has not been used in [2]. This cut improves the ratio of signal to the background 
by about two times. 

The number of events corresponding to signal and background after cuts is presented in the Tables 
2,3 for y/s = 2 TeV and y/s = 4 TeV respectively. One can see that using 3-standard deviation criteria 
one can find Higgs up to mass ~ 100 GeV at y/s = 2 TeV and ~ 120 GeV at y/s = 4 TeV. These 
numbers are obtained if we assume that the efficiency of double b-tagging system e b is about 50% 
for bb pair registration. A rather high e b is very important for the task of the Higgs search, because 
as it was mentioned in [1, 2] there is a huge background from the events with light jets production: 
pp — > Wjj + X; pp — >■ Zjj + X, which has to be suppressed by efficient b-tagging. 

One can see that with the decreasing Higgs boson mass the contribution of the bt process (13) 
to the total background becomes more significant. It should be noted also that Wbb background 
increases approximately by factor two with the increase of y/s from 2 to 4 TeV as well as the WH ( 
or ZH ) signal. At the same time the Zbb cross section becomes about 3-4 times higher because the 
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gg — > Zbb subprocess contribution grows about 6 times. That is why the increasing of the collider 
energy from 2 to 4 TeV extends Higgs mass interval only up to 120 GeV. 

In Fig. 7,8 we demonstrate the effective 66-invariant mass distribution for the cases Mjj = 80 GeV, 
y/s = 2 TeV and Mjj = 120 GeV, yfs = 4 TeV respectively without any cuts and with cuts (14)-(15) 
applied. All event numbers correspond to L = 1000 pb -1 and a clear Higgs peak can bee seen in 
Fig. 7(b) and 8(b). 

However the distributions in Fig. 7,8 and event numbers shown in the Table 2,3 have been obtained 
without real b— quark fragmentation into jets. This question will be considered in the nearest future. 

Conclusions 

The present paper continues the study of the Higgs boson search possibility at the TEVATRON 
energies. In comparison with the results of previous considerations [1,2] the detailed analysis of 
the 4-fermion final state reactions (3) - (5) is presented. Among the complete set of Feynman 
diagrams the main contribution to the final states comes from the Higgs signal WH and ZH, from 
the backgrounds WZ } ZZ and gluonic background Wbb, Zbb with subsequent decays of W, Z bosons. 
It is shown that the important contribution to the background is connected with single top quark 
production subreactions qq — > tb or tb with the top decay to Wb. Such a background has not 
been considered in previous papers [1,2]. As a result of various distributions for Higgs signal and 
background a number of optimal cuts (14), (15) has been found. These cuts suppress background 
for more than by factor 40 while the signal only about by factor 2.5. The final effective bb - pair 
mass distribution allows to extract the signal from the Higgs with mass up to about 100 GeV at 
yfs = 2 TeV and up to about 120 GeV at yfs = 4 TeV for the integrated luminosity L = 1000 pb -1 . 

We should stress also that increase of the collider energy by a factor two extends the Higgs mass 
interval not more than 20%. This fact is connected with the high magnification of the gluonic 
background with energy. Therefore for the task of the Higgs search the most important point is the 
increase of the collider luminosity. 
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Figure Captions 



Fig.l - The total cross section for pp — > (7z/, z/z/, //) + bb + X processes (including the Higgs boson 
signal) versus Higgs boson mass for yfs = 2 and 4 TeV. 

Fig. 2 - Different H° signal event distributions for the both mechanisms (f ), (2) and L = fOOO pb -1 , 
for the case M H = 80 GeV and yfs = 2 TeV. 

Fig. 3 - Different H° signal event distributions for the both mechanisms (f ), (2) and L = fOOO pb -1 , 
for the case M H = 1 20 GeV and yfs = 4 TeV. 

Fig. 4 - the diagrams for QCD 2 — > 3 type main background subprocesses qq — > Wbb } qq — > Zbb 
gg — > Zbb which were included into our generator. Note, that diagrams for another quark flavors 
were also included but are not shown here. 

Fig. 5, 6 — different distributions for background processes (like for the signal) for y/s = 2 TeV 
and yfs = 4 TeV respectively: pp — > WZ — > !vbb, pp — > Wbb — > !vbb, pp — > ZZ — > bbll^v 1 ), 
pp — > Zbb — > bbll^p). 

Fig. 7,8 — the effective 66-invariant mass distribution for the cases Mh = 80 GeV, yfs = 2 TeV 
(Fig. 7) and Mh = 120 GeV, yfs = 4 TeV (Fig. 8) respectively without any cuts ( 7(a), 8(a) ) and 
with cuts (14)-(15) applied ( 7(b), 8(b) ). All events number correspond to the L = 1000 pb -1 . 



Table Captions 

Table 1 — Higgs production rates for processes (1) and (2) with signature (3), (4) and (5) at 
luminosity L = 1000 pb -1 . 

Table 2 — Number of 4-fermion final state events for Higgs boson production (process (1), (2)) 
and background (processes (11)-(13)) for yfs = 2 TeV with signature (3)-(5) and cuts (14), (15) at 
L = 1000 pb -1 and bb pair tagging efficiency 50%. 

Table 3 — Number of 4-fermion final state events for Higgs boson production (process (1), (2)) 
and background (processes (11)-(13)) for yfs = 4 TeV with signature (3)-(5) and cuts (14), (15) at 
L = 1000 pb -1 and bb pair tagging efficiency 50%. 



7 



Table 1: 



M H (GeV) 


pp — >■ W ± H - 


■> 2b- jet + lepton + z/ 


pp — > ZH — > 


2b- jet + 2-lepton 




\j s = Z lev 


\j s = 4 lev 


\j s = Z lev 


\j s = 4 lev 


60 


345 


619 


199 


355 
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307 
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171 


100 


83 


170 
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96 
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Table 2: 



M H (GeV) WE ZH Wbb + W Z + bt(bt) Zbb + ZZ 



60 


35 


27 


53 


+ 1 


+ 8 


41 


+ o 


80 


27 


19 


45 


+ 43 


+ 12 


40 


+ 28 


100 


17 
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+ 44 


+ 12 


32 


+ 29 


120 


9 


7 


26 


+ 1 


+ 13 


25 


+ 1 


140 
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18 


+ o 


+ 11 


17 


+ o 
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Table 3: 



M H (GeV) WH ZH Wbb+WZ + bt(bt) Zbb + ZZ 



60 
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+ 3 


+ 18 


122 


+ 1 
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30 


74 


+ 83 
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+ 69 
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+ 20 


69 
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